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CN ■ ABSTRACT 
> ■ 

oo ; 

I We study the polarization anisotropy in the Cosmic Microwave Background (CMB) 

I resulting from patchy reionization of the IGM by stars in galaxies. It is well known 

' that the polarization of the CMB is very sensitive to the details of reionization, 

■ including the reionization epoch and the density fluctuations in the ionized gas. We 

■ calculate the effects of reionization by combining a semi-analytic model of galaxy 
Q^. formation, which predicts the redshifts and luminosities of the ionizing sources, with 
Q , a high resolution N-body simulation, to predict the spatial distribution of the ionized 
|H I gas. The models predict reionization at redshifts z ~ 5 — 10, with electron scattering 
^ I optical depths due to reionization ~ 0.014 — 0.05. We find that reionization generates 
^ ' a peak in the polarization spectrum with amplitude 0.05 ~ 0.15/iK at large angular 

V~j ■ scales (/ 3). The position of this peak reveals the size of the horizon at reionization, 

j_j ■ whilst its amplitude is a measure of the optical depth to reionization. On small 



scales (/ >6000), reionization produces a second-order polarization signal due to the 
coupling of fluctuations in the free electron density with the quadrupole moment of 
the temperature anisotropy. Careful study reveals that this coupling generates equal 
second-order polarization power spectra for the electric and magnetic modes, with 
amplitude ~ lOnK. This amplitude depends strongly on the total baryon density fib 
and on the spatial correlations of the free electron density, and weakly on the fraction 
/esc of ionizing photons able to escape their source galaxy. The first- and second-order 
signals are therefore sensitive to different details of how the reionization occurred. 
Detection of these signals will place important constraints on the reionization history 
of the Universe. 



Subject headings: cosmology: theory cosmologyxosmic microwave background - 

- galaxies: formation intergalactic medium 
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1. INTRODUCTION 

The secondary anisotropies in the Cosmic Microwave Background (CMB) provide a laboratory 
for the study of the epoch of reionization in the Universe. With the rapid improvement in CMB 
experimental sensitivity and resolution, it is timely to discuss the signals introduced by the process 
of reionization, which are at the /iK level or below. In this paper, we concentrate on the CMB 
polarization in models with inhomogeneous reionization. It is well known that the primordial 
CMB polarization is generated at the recombination epoch through Thomson scattering of the 
quadrupole of the temperature anisotropy. The same mechanism operates at the reionization 
epoch and again distorts the shape of the polarization power spectrum. It has been pointed 
out by Ng &; Ng (1996) and Zaldarriaga (1997) that this leads to the polarization anisotropy 
being suppressed on small scales but enhanced on large scales. The degree of suppression and 
enhancement typically depend on the optical depth to the recombination epoch produced by 
reionization. Although the amplitude of the polarization anisotropy is estimated to be much 
smaller than that in the temperature, it is advantageous to consider the polarization signal since 
it is generated when photons and electrons scatter for the last time. The temperature fluctuations 
are produced both by scattering between photons and electrons, and by the gravitational redshift 
or blue-shift caused by the evolution of the gravitational potential between the last-scattering 



epoch and the present Sachs &: Wolfe 1967. On the other hand, the polarization is affected by 



the gravitational field only through gravitational lensing Zaldarriaga &i Seljak 1998| , which causes 



some smearing of the power spectrum, and also some mixing of the "electric" and "magnetic" 
parity modes of the polarization. In Benson et al. (2001; hereafter Paper I), we presented a 
detailed calculation of inhomogeneous reionization of the intergalactic medium (IGM), using a 
semi-analytic model of galaxy formation to predict the ionizing luminosities of stars in galaxies at 
different redshifts, and coupling this with an N-body simulation to predict the spatial distribution 
of the ionizing sources and of the ionized gas. In this paper, we use the same model to predict the 
effects of reionization on the polarization of the CMB. 

Reionization produces interesting effects on CMB temperature and polarization anisotropies 
at both first and second order in the perturbations. At first order, the effects of reionization are 
the same as for an IGM with spatially uniform density and ionized fraction. Density fluctuations in 
the free electrons around the reionization epoch produce CMB anisotropies only at second order. 



as in the Vishniac effect Vishniac 1987. These second-order anisotropies are small in amplitude, 
but nonetheless dominate over the first-order anisotropies on small angular scales. They are 
thus cosmologically interesting as a probe of structure present at reionization. Historically, work 
on the second-order effects began by considering the temperature anisotropies in the case of a 
homogeneous reionization of the IGM, with the fluctuations in the free electron density being 
assumed to follow the variations in the total matter density, the so-called density modulation model 



Vishniac 1987| , |Jaffe fc Kamionkowski 1998| . However, in any realistic model, the reionization 
is expected to be patchy or inhomogeneous, with some regions already being fully ionized while 
others are still neutral, and the ionized regions growing until they encompass the whole IGM. 
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Shapiro & Giroux (1987) and Sato (1994) investigated the evolution of spherical ionization fronts 
in a homogeneous IGM using analytic methods. More recently, numerical simulations of radiative 
hydrodynamics have also shown how reionization proceeds in an inhomogeneous way Gnedin 2000| . 



Gruzinov & Hu (1998) investigated second-order temperature anisotropies resulting from patchy 
reionization using a simple analytical model in which the universe is reionized in uncorrelated 
spherical patches. More realistic calculations which relax the assumption of uncorrelated spherical 



ionized regions have since been carried out by various authors, using either full simulations [Gnedin 



fc Jaffe 2000| or semi-analytic methods (Paper I, Springel, White, & Hernquist 2000, Bruscoli et al. 
2000, Valageas, Balbi, & Silk 2000). So far, these works have concluded that second-order effects 
make significant contributions to the temperature anisotropy on small scales (/ >3000), with peak 
amplitude in the range 0.1//K to 1/iK. 

While the second-order temperature anisotropies are well studied, those of the polarization 
have received relatively little attention. Seshadri & Subramanian (1998) and Weller (1999) 
discussed this effect in density modulation and patchy reionization models, respectively. Here, 
we consider second-order polarization anisotropies in a more realistic model than has previously 
been possible. In our model, reionization results from photoionization by stars in galaxies, and the 
spatial fluctuations in the free electron density are the combined effect of fluctuations in the total 
density and fluctuations in the ionized fraction, thus combining the "density modulation" and 
"patchy reionization" effects. We calculate these effects in the same way as in our previous work 
(Paper I), by combining a semi-analytic model of galaxy formation with an N-body simulation of 
the distribution of dark matter in the universe to determine the distribution of ionized regions. 
In this paper, we consider only the scalar mode for the primordial fluctuations. In the case of 
flrst-order perturbation theory, the scalar mode fluctuations produce only the electric (E-) mode 
of polarization Zaldarriaga &: Seljak 199^ . For the second-order perturbations, however, this 



may not be the case. Since the polarization anisotropies are produced by the coupling between 
the large-scale primordial CMB temperature quadrupole anisotropies and the small-scale density 
fluctuations in the ionized medium, axisymmetry is broken and the magnetic (B-) mode is also 
produced Hu 2000| . In addition to the primordial CMB temperature quadrupole anisotropy, there 



is also a kinetic temperature quadrupole in the rest-frame of the scatterers, produced by the 



quadratic Doppler effect from the peculiar velocity of the ionized medium Sunyaev &: Zel'dovich 



1980 . Coupling between this kinetic quadrupole and the electron density fluctuations can likewise 
produce a second-order contribution to the polarization, but this contribution is much smaller 
than the previous one Hu 200C| , so we neglect it in the present paper. 



This paper is organized as follows. In §2, we present the Boltzmann equation for the first- 
and second-order polarization anisotropies. We then derive the E- and B-mode power spectra by 
use of the Boltzmann equation. In §3, we give a brief overview of our model for inhomogeneous 
reionization and show our numerical results for both the first- and second-order polarization power 
spectra. §4 is devoted to our conclusions. Throughout this paper we work in units where c = 1. 
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2. SECOND ORDER EFFECTS FROM THE REIONIZATION EPOCH 

2.1. Boltzmann Equation and Second Order Polarization 

The evolution of the temperature perturbation, A7^(x, n, r), and the polarization perturbation, 
Ap(x, n, r), are governed by the Boltzmann equation |Bond &z Efstathiou 1984 , , . Here x is the 



comoving coordinate, r = J dt/a the conformal time, where a is the expansion factor normalized 
to unity today, and n the direction of photon propagation. Following Zaldarriaga & Seljak 
(1997), to calculate the polarization perturbation Ap we work in terms of the perturbed Stokes 
parameters AQ±iij (see §2.2). In this paper, we do not include tensor perturbations in the 
metric, and we derive the Boltzmann equation in the conformal Newtonian gauge. Readers 
interested in the synchronous gauge or the transformation between these two gauges are referred 
to Ma & Bertschinger (1995). In the Newtonian gauge, the perturbations are specified by two 
scalar potentials, (p ^-nd ip, which play the role of the gravitational potential and the fractional 
perturbation to the spatial curvature, respectively. The Boltzmann equation then states that 

At + = 4>- niditp + n^arair) -At + A^o + n^vu " ^ E >1"(n)n(") , (1) 



-AQ±it/ + yy E ±2>r(n)n("^M , (2) 

where sY[^ is a spherical harmonic with spin-weight s, and 11^™-* is defined in terms of the 
quadrupole components of the temperature and polarization perturbations as 

n(") (x, r) EE aJ^) (x, r) + 12 V6A^ (x, r) + uVqa'lI (x, r) . (3) 

Here the derivatives are taken with respect to the conformal time, is the velocity of baryons, ne 
is the free electron number density, and o"t is the Thomson cross section. We have also expanded 
the perturbations in the radiation field in spherical harmonics with appropriate spin weight s 
(see also §2.2) 



AT(x,n,r) = E(-^)V4^(2Z + l)A?;^)(x,r) Yr{ii) 

Im 



AQ±,f;(x,n,r) = Y^{-zy^A7r{2l + l)Ag)(x,r) ±2Yr{n). (4) 

Im 



When we calculate the perturbations for a mode with wavenumber k, we define the sl^"^'s in 
a coordinate system for which the z-axis is parallel to k. To calculate the first order effect for 
scalar modes, we can set m = 0, as in Zaldarriaga & Seljak (1997), due to the axisymmetry of 
the radiation field around the mode axis for this case. However, their expansion is not valid for 
studying the second order effect, in which the rotational symmetry around the wavevector is 
broken by coupling to other modes. It is very important to take into account m ^ modes for the 
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calculating the second order effect, otherwise the (artificially imposed) axisymmetry guarantees 
that no magnetic mode can be generated. This is why Weller (1999), who assumed axial symmetry 
following Seshadri &: Subramanian (1998), obtained only the electric mode of polarization for the 
second order effect. As was found by Hu (2000), who employed the weak-coupling approximation, 
the magnetic mode of polarization is also generated in the case of the second order effect. Thus, 
the expansion in equation (^) is more general and useful for our calculation. 

On small scales, as pointed out by Kaiser (1984), the contribution to secondary temperature 
anisotropics from homogeneous reionization tends to cancel along the line of sight. Thus, the 
first order effect on the temperature and polarization anisotropics from the reionization epoch is 
suppressed at small angles. Hereafter, we concentrate on polarization and develop an equation for 
the visibility- modulated effect which is the dominant source on small scales. 

Firstly, we write the inhomogeneous distribution of the free electron number density as 

ne(x,r) = ne(r)[l + 5e(x,r)], (5) 

where 6e is the fluctuation in the electron number density and "denotes the background quantity. 
Then, equation (|2|) can be rewritten in terms of Fourier modes as follows, 

^Q±iu{K t) + iA;/iAQ±i(/(k, r) = neCrTa{T) ( -AQ±iu{k, r) - R±{k, r) 

^ ^ (6) 



+ E\/t ±2l'2'"(fi)(n('")(k,r) + 5('")(k,r))) , 

m / 



where 5('")(k,r) and R±{k, r) are the Fourier modes of the coupling of (5e(x, r) to n('"^(x, r) and 
^Q±iC/(x, r), respectively. In other words, 

5™(k,r) = J dVc(k-p,r)n(-)(p,r), (7) 

R±{K t) = j d3p<5e(k - p, r)AQ±it;(p, r). (8) 
The polarization perturbations at the present time can be obtained by integrating the Boltzmann 



equation (g) along the line of sight Zaldarriaga &: Seljak 1997 
AQ±.f/(k,/z,ro) = / dre-*'=(^«-^)^g(r) 



X 



5 

where = k • n, and g{T) is the visibility function defined as 



±2ir(fi)K-)(k,r) + 5(™)(k,r)l -i?±(k,r)), (9) 



5(r)^-^e'^(^")-«W, (10) 
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with k{t) = J^" dr'arieaT being the electron-scattering optical depth. The visibility function 
has a simple physical meaning, being the probability that a photon had its last scattering at 
epoch r and reaches the observer at the present time, tq. Equation can then be given a 
simple interpretation, since n('")(k,r) acts as a first order source term while i?±(k) and 5('")(k) 
are the contributions from the second order effect. To simplify the calculation, we neglect 
i?±(k) and A^^l in 5('")(k) because the temperature perturbations typically dominate over the 

polarization perturbations. Furthermore, we assume the first order temperature quadrupole Aj^^ 
is uncorrelated with 5^ due to the fact that the dominant contributions to these come from large 
and small scales, respectively. That is, we regard the source term for the second order polarization 
as 

5(-)(k,r) ^ <^e(k,r) j d3pA5^)(p,r) 

^ <5e(k,r)Qf)(r), (11) 

with g^"V) defined as = /d^P^S^^ (p, r) for convenience. Recall that the scalar mode 

in linear theory only generates the m = component in the p-basis, i.e. AP^Y^{(5,a), where (3 
and a are the polar and azimuthal angles defining n in this basis. Using the addition theorem, we 
can project the component in the p-basis onto the k- basis (see Ng & Liu 1999), 



^il™*(n)yr(p) = J^i^AA"), (12) 

It then follows that 



^^"^ = Vy / d3pA?^(p)yr(p)- (13) 

Finally, the solution for /S.Q±iu becomes 



AQ±,[;(k,n,ro) = J^ rdTe^'=(-°--)^5(r)E ±2!^ (fi)^^"^^ (k, r), (14) 

where X('")(k,T) equals n(o)(k, r) or 5*-™^ (k, r) for the first and second order contributions 
respectively. 



2.2. Stokes Parameters and the Power Spectrum 

We follow the approach of Zaldarriaga & Seljak (1997) of calculating the polarization in terms 
of the electric and magnetic modes. We start with the description of the polarization perturbation 
in terms of the Stokes parameters Q and U. If we consider a wave traveling in the z direction, 
Q is the difference in intensity polarized in the y and x directions, while U is the difference in 
the (x -|- y)/\/2 and (x — y)j\pl directions. The circular polarization parameter, 1/, cannot be 
produced by scattering, so we will not mention it further. Polarization is more complicated than 
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temperature because the values for Q and U depend on the choice of coordinate system. Under a 
right-handed rotation through an angle around the z axis the temperature is invariant, while Q 
and U transform according to 

{Q±iU)'{n) = e^'''^{Q±iU){h). (15) 

Thus the functions Q ziziU have spin ±2, and should be expanded in spherical harmonics of spin 
±2. 

It is possible to produce a rotationally invariant measure of the polarization field if we 
introduce the spin raising and lowering operators 3 and d Newman &: Penrose 196(^ 



= -(sin6')'' —+icsc9— (sm9) "rj, 
[ou ocpj 

&n = -(sin^)-" -7^ -i CSC 9— {sm9yi], (16) 
lo9 ocpj 

with ?7 as a spin-s field. If they act on the spin-s spherical harmonics, we have 

~dsYlm = -[{l + sKl-S + l)]K^lYlm. (17) 

The new bases of the rotationally invariant polarization field, called the electric mode and 
magnetic mode, are defined as 

As = -'-(6^AQ+iu-&^^Q^^u). (18) 



2 

The values of A^; and A^ at a particular direction in the sky are independent of the coordinate 
system used to define them (unlike Aq and Ajj)- We work in the coordinate system where k \\ z 
and define Q > (Q < 0) in the direction eg{e^). For a scalar mode in first-order perturbation 
theory, the radiation field is axisymmetric around the wavevector. The polarization is produced by 
scattering of the quadrupolar component of the radiation field, and so has only a Aq component, 
but no Ajj component (which would correspond to a polarization angle at 45° to the eg - 
axes), thus Aq^hj = Aq^hj = Aq. Furthermore, Aq has no (/^-dependence, so 3 ^ = 6 ^ (eq. 
16). As a consequence, A^; = —d ^Aq = —6 ^Aq and A^ = 0. This result is important because 



it shows that no magnetic mode can be produced by scalar modes (density perturbations) in 
first-order perturbation theory. However, as we will show later, the same situation does not occur 
in the second-order contribution to the polarization because the coupling between modes breaks 
axisymmetry. 

Again, following Zaldarriaga k. Seljak (1997), the power spectra of the electric mode and 
magnetic mode can be defined as 

^(^'^)' " ^rrijp^^ I < = ro)p >, (19) 
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where a|^^^^^ can be extracted using 



(20) 



To calculate equation (|20D, we can apply the spin raising and lowering operators (eq. |16| or 
17 ) to the polarization perturbations in equation (^) twice. But the plane wave e^K'ro-'r)p. j^ggjf 
carries an angular dependence, thus we expand the plane wave in a series of spherical waves using 
the Rayleigh equation 



5:(-^)74vr(2/ + l)j,[A;(ro-r)]yAA) 



(21) 



where ji[k{To — r)] is the spherical Bessel function. Furthermore, we calculate the product of two 



spherical harmonics with spin by using the Clebsch-Gordan relation Sakurai 1985 



si^li h ~ / , < ii,i2,mi,m2\i,m > 



Ims 



X < li,l2, —Si, —S2\l, —S > 



I Air 
21 + 1 



s^l ■ 



(22) 



The expression for A|^''^^^(k, r) then becomes 

3 ^(/ + 2)' 



A 



-^/2lTl rdr5(T)X(")(k,r)rg)^),(fcr 
J 



(23) 



with r = tq — T. In Table 1, we list the results of T^^^^^{kr), in which the orthonormality relation 

^U'Smm'i (24) 

and the recursion relations of spherical Bessel functions 



1 



21 + 



j[ji-i{x) + ji+i{x)] 



(25) 



have been used. Finally, the power spectrum can be obtained by substituting equation (|23| ) into 



C^E,B)l = (4vr)^^|±|jE / k'dkllj dr,(r)xM(k,r)Tg;),,(^r) 



(26) 



In the case of the second-order effect, eq.(|26|) can be rewritten as an expression involving the 
power spectrum of the fluctuations in the free electron density (5e(k, r)(5*(k, r')), connecting the 
two different times r and r'. When we evaluate eq.(^) numerically, we neglect the effect of the 
difference in times in the power spectrum, and so replace it by {6e0^,f)6*{]<.,f)) , evaluated at a 
single time f which is an average of r and r'. This should be an adequate approximation if the 
dominant wavelengths in eq.(26) are small compared to the Hubble radius. 



m 



'El 



±1 



■21+1 Y 24(/+2)! V 



-I- / (-0' 



±2 



(f+2)(;+i) 

2«-l 



2i+3 



6 



(2;+i)(;-i)(;+2) 

(2/-l)(2«+3) 



zt) - (/ + 2)(/ + + - 



■(/-l)ii+i(A;r)] 



Table 1: T^^^^^y in equation (|2 



3. NUMERICAL RESULTS 



To calculate the polarization anisotropy spectrum using the results discussed above, we use 



the publicly available code CMBFast peljak fc Zaldarriaga 1996| . We modify the ionization history 
in this code to follow the more realistic case from our previous work (Paper I), in which the 
reionization history of the universe is determined by a semi-analytic model of galaxy formation. 
The semi-analytic model is that of Cole et al. (2000), which includes the following processes: 
formatinon and merging of dark matter halos through hierarchical clustering; shock-heating and 
radiative cooling of gas within these halos; collapse of cooled gas to form galactic disks; star 
formation in disks and feedback from supernova explosions; galaxy mergers; chemical evolution of 
the stars and gas; and luminosity evolution of stellar populations based on stellar evolution codes 
and model stellar atmospheres. The model has been shown by Cole et al. to agree well with a 
wide range of observed properties of galaxies in the local universe. In Paper I we used this model 
to calculate the ionizing luminosities of galaxies at different redshifts, including the effects of 
absorption by interstellar gas and dust on the fraction of ionizing photons escaping, and followed 
the propagation of the ionization fronts around each galaxy. To find the ionizing luminosity, we 
first calculate the rate at which ionizing photons are being produced by stars in the galaxy, then 
apply an attenuation due to dust, and finally allow a fraction /esc of the remaining photons to 
escape into the IGM. The mass of ionized hydrogen in each spherical ionization front is found by 
integrating the equation ^hapiro &: Giroux 1987 



1 dM 
mil dt 



(2) —'^ 

S{t)-ali'a /clump 



M 
mn' 



(27) 



where nn is the comoving mean number density of hydrogen atoms (total, Hi and Hll) in the 
IGM, nin is the mass of a hydrogen atom, a{t) is the scale factor of the universe normalized to 
unity at z = 0, i is time (related to the conformal time by dt = dr/a), S{t) is the rate at which 

(2) 

ionizing photons are being emitted and a^j is the case B recombination coefficient. The clumping 
factor /clump = {Pigm)/p'igm gives the effect of clumping on the recombination rate of hydrogen in 
the IGM. A larger /dump increases the recombination rate resulting in a delay of the reionization 
epoch. We use the clumping factor f^i^^^ as defined in Paper I. By summing over the ionized 
volumes due to all galaxies, we can calculate the fraction of the IGM which has been reionized at 



- 10 - 



any redshift. Having this reionization history, we can then obtain the first order power spectrum 
of the polarization anisotropics without making any assumptions about the spatial distribution of 
the ionized gas. We consider only the scalar mode of the primordial perturbations, for which the 
radiation field for each k-mode is axially symmetric. Thus, we set m = in equation (26) and find 
no magnetic mode can be produced (see Table. 1). 

The second order effect is more complicated, as the source term contains 6c and i.e., 
the m = 0, ±1 and ±2 components of the temperature quadrupole must be considered. The time 
evolution of the temperature quadrupole components with different m in the k-basis is obtained 
using equation (^), and A^?] is calculated using the modified CMBFast code. This leaves only the 
power spectrum of the electron density fluctuations 6e unknown, which should be calculated from 
the distribution of ionized regions with different sizes and shapes and with a correlated spatial 
distribution. An exact calculation would require a high resolution numerical simulation with gas 
dynamics and radiative transfer included, and so would require very large amounts of computing 
time. Here we instead calculate the density field 5c of the ionized gas using the simpler approach 
of Paper I, in which the semi-analytic galaxy formation model is combined with a high-resolution 
N-body simulation of the dark matter. The simulation volume, which is a box of comoving length 
141.3/i~^Mpc and contains 256^ dark matter particles each of mass 1.4 x 10^^h~^MQ, is divided 
into 256^ cubical cells of comoving size 0.55/i~^Mpc. Each dark matter halo in the simulation 
containing more than 10 particles is populated with galaxies according to the semi-analytic model, 
and the sum of their ionizing luminosities placed in a source at the center of mass of the halo. 
Ionizing photons which originate from galaxies in lower mass halos (which are not resolved in the 
simulation) are added in by assuming that these unresolved halos trace the remaining mass of the 
simulation (i.e. that mass which is not part of a resolved halo). In Paper I we demonstrated that 
the locations of unresolved sources do not significantly affect the resulting anisotropy spectrum. 
We determine which regions of the simulation box become ionized by using one of the five toy 
models A-E listed below, which span the likely range of possible behaviour. In all cases, the total 
mass of ionized hydrogen in the simulation box is the same, and is forced to equal that for a 
large-scale homogeneous distribution with the specified IGM clumping factor /dump) calculated by 
use of equation (^). We assume that the total gas density in the IGM traces the dark matter, 
thus neglecting the effects of pressure in the IGM. Hu (2000) has shown that this is a good 
approximation for anisotropies with I ^10^. 

Model A (Growing front model) Ionize a spherical volume around each halo with a 
radius equal to the ionization front radius for that halo calculated assuming a large-scale uniform 
IGM. Since in the simulation the IGM is not uniform, but is assumed to trace the dark matter, 
and also because some spheres will overlap, the ionized volume calculated in this way will not 
contain the correct total ionized mass. We therefore scale the radius of each sphere by a constant 
factor, /, and repeat the procedure. This process is repeated, with a new value of / each time, 
until the correct total mass of hydrogen has been ionized. 



Model B (High density model) In this model, we ignore the positions of halos in the 
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simulation. Instead, we simply rank the cells in the simulation volume by their density. We then 
completely ionize the gas in the densest cell. If this has not ionized enough Hi then we ionize the 
second densest cell. This process is repeated until the correct total mass of Hi has been ionized. 

Model C (Low density model) This is like model B, except that we begin by ionizing the 
least dense cell, and work oTir way up to cells of greater and greater density. This model mimics 
that of Miralda-Escudc, Hachnelt & Rees (2000). 

Model D (Random spheres model) As Model A, except that the ionized spheres are 
placed at completely random positions in the simulation volume, rather than on the dark matter 
halos to which they belong. By comparing to Model A this model allows us to estimate the 
importance of the spatial clustering of dark matter halos. 

Model E (Boundary model) Ionize a spherical region around each halo with a radius equal 
to the ionization front radius for that halo. This may ionize too much or not enough Hi depending 
on the density of gas around each source. We therefore begin adding or removing cells at random 
from the boundaries of the already ionized regions until the required mass of Hi is ionized. 

Prom the results of the simulation, we calculate Se as 

Xe(x, t){1 + S) 

4(x,r) = 1, (28) 

where Xe is the ionized fraction (which we take to equal 1 in ionized regions and in neutral 
regions) and 6 is the dark matter overdensity (i.e. we assume that fluctuations in the gas density 
follow those in the dark matter). Here Xe is defined as the mass-averaged ionized fraction in the 
IGM. 

We fix the cosmological parameters to be = 0.3, Aq = 0.7, Hubble constant Hq = 70 
km/s/Mpc and as = 0.9. We will consider the effects of varying /esc and on the polarization 
anisotropics and also examine the polarization power spectrum in all five toy models for the 
distribution of ionized regions. In Fig. 1 we plot the visibility functions of the ionization histories 
in our simulation for different fesc (panel (a)) and different 0,^ (panel (b)). We find that the 
visibility function depends strongly on /esc, and also on Q,^. The value of J^b determines the 
cooling rate (and so star formation rate) in our model of galaxy formation and also alters the 
recombination rate in the IGM, and so affects the time at which reionization occurs (panel (b)). 
Note that when we vary Jib and /esc, we also vary the fraction of brown dwarfs in the IMF 
used in the galaxy formation model, in such a way that the model always fits the "knee" of the 
observed Ha luminosity function of galaxies at ^ = (see Paper I for more details). Therefore the 
production of ionizing photons does not simply scale with 

In Fig. 2, we show how the first-order polarization anisotropics are affected by scattering 
by free electrons at the reionization epoch, calculated using the modified CMBFast code. We 
find, first, that rescattering at the reionization epoch generates a new anisotropy at large scales 
because the horizon has grown to a much larger size by reionization than it had at recombination 
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(z ~ 1100). More specifically, the location of the new peak reveals the horizon size at last 
scattering, and its height reveals the duration of last scattering, i.e., this new peak is sensitive 
to the optical depth produced by reionization. To see how the distortion of the primordial 
polarization depends on the optical depth Kj back to the start of reionization, we plot the resulting 
power spectrum for various values of /esc in panel (a). We find that the boost in the large-scale 
power depends strongly on the value of /esc- (The optical depth is 0.034, 0.017 and 0.014 for 
/esc = 1.0,0.1 and 0.05 respectively, for fib = 0.02.) The location of the new peak depends on the 
reionization epoch, which is characterized by the redshift Zi of the peak of the visibility function 
as shown in Fig. 1. From numerical simulations of the Boltzmann code for different reionization 
epochs and cosmological parameters, we obtain a fitting formula for the peak location ^peak as a 
function of the reionization epoch Zji 



This fit is consistent with the peak locations shown in Fig. 2. On the other hand, the first-order 
polarization fluctuations are suppressed on small scales by rescattering since a fraction of the 
photons coming from the recombination epoch is scattered after reionization. The tiny optical 
depth for rescattering causes little erasure of power on small scales, but the new peak reaches an 
amplitude of ~0.05~0.1/iK. From Fig. 2(b), we can see that the suppression also depends on Kj. 

Let us discuss now the second-order effect. In Fig. 3, we show the contribution to Cei and 
Cbi from each value of m for Model A. For m = 0, only the E-mode is generated because = 
(see Table. 1). Our numerical results show that the E-mode (the B-mode) from m = ±2 (m = ±1) 
dominates the one from ±1 (±2) by more than four orders of magnitude. However, the total 
power spectra of the E- and B-modes are almost exactly the same. The difference is less than 
10~^, which may be caused by numerical errors. These results are consistent with previous work 
by Hu (2000), in which the weak-coupling approximation was employed. (The weak-coupling 
approximation involves making an analytical approximation for the line-of-sight integrals, while in 
eq.(^6|) we instead calculate the integrals over r and k numerically.) The reason for the equality 
of the total E- and B-mode power spectra is essentially that the first-order quadrupole whose 
scattering produces the polarization is dominated by large scales, and so has a random orientation 
relative to the small-scale fluctuations in the electron density (c.f eqn.(|ll|)). Scattering of the 
quadrupole by the small-scale fluctuations therefore on average excites E- and B- modes equally. 

Numerically, the amplitude of the second-order signal is found to be about 10~^K. This 



can also be understood as follows: from (eq. the amplitude of the second order polarization 
anisotropics is roughly A.Q±iif ~ / dTg{T)Q2Ue ~ nQ^^^'de- In our calculation k, and 5e 
are typically on the order of 10~^, 10~^K and 10~^, respectively. The order-of-magnitude of the 
second-order contribution then follows. 




(29) 




The relation between the angular wavenumber / and the comoving wavenumber k of the 
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density fluctuations at a given redshift is 

I ~ kr{z) = k / , ^ ^, ■ (30) 

where r{z) is the standard radial coordinate distance in the Robertson- Walker metric, and the 
final expression assumes a flat cosmology. For our standard choice Qq = 0.3, Aq = 0.7, we 
find r{z) « 6000/i~^Mpc for reionization at ~ 5 — 10. Thus, electron density fluctuations at 
reionization on comoving scale Ax produce second-order anisotropies at I ~ 6000 {h~^Mpc/Ax). 

In Fig. 4, we plot the second-order power spectrum of the polarization in the five toy models 
with fixed extreme escape fraction /esc = 1 and Uf, = 0.02. Although the shapes of the curves are 
all very similar, their amplitudes are different. Note that the reduction in power above I ~ 10, 000 
is artificial and due to the limited resolution of the N-body simulation we use (the density field 
of the ionized gas is calculated on a grid with cell size 0.55/i~^Mpc, corresponding to Z ~ 10^). 
On the other hand, the finite size of the simulation box (256/i~^Mpc) affects the power spectrum 
for I below a few hundred. We see that the amplitude of the power spectrum around the peak 
{I 10, 000) varies by a factor Ri 2.5, depending on which of the models A-E is used. In Paper I we 
argued that the amplitude of the curves is affected by the strength of the correlations of Sg present 
in each model. As a result, the "high density" model (B) is the most strongly correlated and has 
the highest amplitude, and conversely the "low density" model (C) has the lowest amplitude. 

In Fig. 4, we also compare our results to the analytical toy model of Gruzinov &: Hu (1998), in 
which the reionization is described by three free parameters. In their model, each luminous source 
is assumed to ionize a spherical region with fixed comoving radius R, the first source appears at 
redshift Zi, and new sources turn on at a constant rate until reionization is complete after an 
interval Sz. An artificial assumption is made that luminous sources appear at random locations 
in space, so there are no correlations between the positions of the ionized spheres. Assuming 
that the spheres remain ionized forever, the fractional ionization increases with increasing number 
density of ionized spheres during 5z until the universe is completely ionized. In Paper I, we chose 
R = 0.85/i^^ Mpc, Zi=ll and 6z = 5 in the Gruzinov & Hu model to match the peak in the 
power spectrum of secondary temperature anisotropies predicted for our model E. We choose the 
same parameters here and find that the toy model likewise matches the peak in the spectrum of 
polarization anisotropies for model E. For small I, little power is generated in the Gruzinov &; Hu 
model because, by design, the patches are uncorrelated. 

To further clarify what determines the shape of the second-order anisotropy spectrum in 
our models, we carried out the following additional tests, shown in Fig. 5. The first test was to 
force the ionized fraction Xg to be uniform and equal to the same mean value as before, so that 
fluctuations 5e in the free electron density are then simply equal to fluctuations S in the total 
density. In this case, the angular power spectrum has an almost identical shape to models A-E, 
with amplitude about equal to that for models C, D or E, but lower than that for models A or B 
(by factors of about 1.7 and 2.5 respectively). In the next two tests (labelled "random" in Fig. 5), 
the total gas density was forced to be uniform (i.e. we set S = 0), and put bubbles down at random 
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positions, so that fluctuations in 5e resulted only from the patchiness of the reionization. In one 
case, the bubble radii were chosen from the size distribution predicted by our galaxy formation 
model (this case is thus similar to model D, except that model D had 5 7^ 0). In the other 
case, all bubbles were given the same comoving radius of 0.62/i~^Mpc, which corresponds to the 
mean bubble radius (weighted by bubble volume) predicted by the galaxy model, at the redshift 
corresponding to the peak of the visibility curve. Both of these cases give power spectra with 
shapes (at large scales) similar to the analytical Gruzinov & Hu model, and completely different 
from when fluctuations in 5 are included. In the final test (labelled "clustered" in Fig. 5), we 
again forced the bubble radii to be equal at a given redshift, but placed them on random halos, 
and included the fluctuations in the total gas density. The starting value for the radii in this last 
case was again 0.62/i~^Mpc, but the spheres were then grown by a uniform factor at each redshift 
to produce the correct mean ionized fraction, as in model A. The power spectrum in this case is 
almost the same as in model A, showing that the distribution in bubble sizes in the latter case 
does not have much effect. 

From these tests, we conclude that in our models A-E, the shape of the power spectrum 
on scales large compared to the typical size of the ionized bubbles is determined primarily by 
the correlations in total gas density. However, in the case of patchy reionization, the amplitude 
depends on the spatial distribution of these patches, which produces biasing for the correlations 
in the ionized gas density relative to those in the total gas density, which in turn boosts the 
amplitude of the polarization fluctuations. For our models A-E, the boost in amplitude of the 
power spectrum, relative to the uniform Xe case, varies from about a factor 1 to 2.5, the largest 
boost resulting from the case where the densest cells are ionized first (model B). 

The amplitude of the second-order polarization power spectrum also depends on the value 
of f^b, which affects the visibility function in equation (p3[). In our calculation, fib is set to be 
0.02, which is lower than the current estimates based on Big Bang Nucleosynthesis which imply 
Ob/i^ = 0.020 ± 0.002 [Buries et al. 2001| . In Fig. 6, we increase fib to 0.04 (which results in m 
increasing to 0.048). If the evolution of the ionized regions were the same in both models, we 
would expect the amplitude to increase a factor of four (eq. [l^ . However, the delay of reionization 
resulting from this increase in fibi as was shown in Fig. 1, reduces the factor to 2.8. 

So far we have set /esc = 1-0 only. Whilst this may in fact be a plausible value for high-redshift 
galaxies based on recent observations (Steidel Pettini & Adelberger 2000), it is instructive to 
examine the effects of changing /esc- Using Model A we examine the effects on the second order 
polarization anisotropics of varying /esc- The results are shown in Fig. 7. In contrast to the case 
of first-order anisotropics, the amplitude of the second-order anisotropics increases only weakly 
with increasing /esc, with the extreme value fesc = 1.0 differing from the others by only a factor 
of about 1.6. This is because the source term for the first-order polarization anisotropics is the 
visibility function multiplied by the temperature quadrupole, g{T)/S.T2- Since the temperature 
quadrupole changes slowly with time, we can say that the power spectrum is mainly determined by 
the visibility function. However, the source term for the second order anisotropics is the visibility 
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function multiplied by Se, whose time evolution is roughly proportional to 1/(1 + ^;). Thus, the 
signal in the second-order anisotropics is more strongly weighted to low z, where the differences 
between models with different /esc are less significant. 

We also plot the first-order power spectrum as heavy solid lines in Figs. 4-7. This shows that 
the second order signal dominates over the first order signal on small scales (/ >6000). This will 
be useful for a high resolution experiment since we will be able to distinguish the secondary signal 
from the primordial one. 

4. CONCLUSIONS 

We have investigated the polarization anisotropies created during the reionization epoch, for 
a realistic model of reionization by galaxies, in particular focusing on the second-order effect. It 
is found that the boost in the large-scale polarization anisotropy from the first-order effect is 
very sensitive to the optical depth due to reionization, with larger optical depths giving larger 
boosts. Using a semi-analytic model of galaxy formation to calculate reionization (Paper I), the 
optical depth from the present to the recombination epoch, m, is predicted to be in the range 
0.014 — 0.048, for ilf, in the range 0.02 — 0.04, and the escape fraction of ionizing photons from 
galaxies /esc in the range 0.05 — 1.0. This first-order effect causes a new peak in the polarization 
at large scales with amplitude around 0.05 — O.lb/jK. The position of the peak is determined by 
the size of the horizon at reionization. If this peak can be detected by future experiments such as 
Planck, it would offer valuable evidence to guide our understanding of reionization of the universe. 

We have also studied the second-order polarization effects resulting from the coupling of 
fluctuations in the free electron number density, with the quadrupole of the temperature 
anisotropies. We obtained (^e(x, by combining the semi-analytic model of galaxy formation 
with a high resolution N-body simulation, as in our previous work (Paper I). The semi-analytic 
model determines the average ionized fraction at each redshift, but the spatial distribution of 
ionized gas is based on using the simulation to determine the locations of the ionizing sources and 
density field of the IGM. We determine which regions of the IGM are ionized using one of five 
toy models. We summarize our results for the second-order polarization anisotropies as follows: 
(1) The second order effect dominates over the first-order effect on small scales {I > 6000). (2) 
The B-mode of polarization is induced as well as the E-mode. The angular power spectra of 
these two modes Cei and Cbi are the same. (3) The shapes of the C(^e,b)i are very similar in 
all of the toy models we considered for the spatial distribution of the ionized regions, but differ 
considerably from the toy model of Gruzinov & Hu (1998), having much more power on larger 
angular scales. The reason for this difference is the spatial correlations of the total gas density 
in our model, which produce corresponding large-scale correlations in the density of ionized gas 
also. The shape of C(^e,b)i on large scales is determined mostly by the correlations in the total gas 
density, but its amplitude is sensitive to the geometry of the ionized regions, which determines 
the biasing factor of the correlations in ionized density relative to those of the total gas density. 
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(4) The difference in amplitude between our different toy models for the spatial distribution of 
ionization is large in the second order effect, for example, Models B and C differ by a factor of 2.5. 
Thus, they provide a very important constraint on galaxy formation and reionization, in spite of 
the fact that the shapes of all curves are similar. The cosmological parameters which also affect 
the amplitude of the second-order effect will be determined by forthcoming precise measurements 
of CMB anisotropics by MAP and PLANCK. Therefore we expect that when the second order 
polarization is observed, its amplitude will provide an important signature of the galaxy formation 
and reionization processes. (5) The amplitude also depends on the value of fib- If the evolution 
of the ionized regions were the same for any value of we would expect the amplitude to be 
proportional to 0^. However, the scaling in our model is somewhat weaker than this, because 
increasing the value of Ob delays reionization. In our calculation, the second-order power spectrum 
for an Jib = 0.04 model has an amplitude 2.8 times greater than for an fib = 0.02 model. (6) The 
amplitude of C(^e,b)i depends only weakly on /esc (and so on the redshift of reionization). Since 
the amplitude of the power spectrum for the extreme value of /esc = 1-0 differs from that with 
/esc = 0.1 only by a factor of about 1.5, the results in Fig. 4 and Fig. 6 should be only weakly 
dependent on the true value of /esc- 

While second-order effects can make a significant contribution to the temperature anisotropy 
(~0.1-1//K), the second-order polarization anisotropics generated by inhomogeneous reionization 
are, as expected, very small. We conclude that the signals are typically of order lOnK on arcminute 
scales, comparable to estimates from other authors studying Vishniac-type polarization [Seshadri 



fc Subramanian 1998| or toy models of patchy reionization Weller 1999 . A signal of this amplitude 



is below the detectability limits of the Planck Surveyor mission, which is the most accurate 
experiment in the near future. However, it is worthwhile to note that this second order effect 
generates the same power in the E- and B-modes. Thus a detection of B-mode polarization 
on arcminute scales or below is essential to distinguish the second-order effect from other 
possible effects. Since it provides more information about galaxy formation and the process of 
reionization than the first-order effect, detection of this signal should be a key aim of a post-Planck 
experiment with increased sensitivity and resolution in the next decade. Among currently 
planned instruments, the most useful one for detecting the second-order signal is probably the 
mm-wavelength interferometer ALMA, which will be able to measure CMB fluctuations on 
arcminute scales. Since the second-order signal is very small, polarized Galactic foregrounds, 
from dust (e.g. Prunet et al 1998| .) and synchrotron emission (e.g. Baccigalupi et al 2001 ), 



are a concern. However, the amplitude of these foregrounds on arcminute scales is currently 
unknown, as observational estimates have been obtained only for larger scales [l <100 — 1000). 
The detectability of the CMB polarization signal can thus only be decided by future observational 
work. 
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Fig. 1. — Visibility function for various escape fractions /esc (panel (a)) and baryon fractions Qh 
(panel (b)) in ACDM models with = 0.3, h = 0.7 and A = 0.7. In panel (a), Ob is fixed at 0.02, 
and in panel (b), /esc is set to 1.0. 
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Fig. 2. — The modification of tlie first-order polarization power spectrum by reionization. (a) The 
angular power spectrum of the first-order polarization produced by the reionization histories with 
Ob = 0.02 and different escape fractions /esc- (b) The fractional change in these power spectra 
relative to the model with no reionization. The boost on large scales and supression on small scales 
are due to reionization. The values of /esc = 0.05,0.1 and 1.0 correspond to optical depths due 
to reionization of Ki = 0.014,0.017 and 0.034 respectively, for fib = 0-02. In panel (a), the power 
spectrum is multiplied by the square of CMB temperature Tcmb = 2.726iir. 
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Fig. 3. — Contributions to the angular power spectrum of polarization from the second-order effect 
from different values of m. The curves shown are computed using Model A with /esc = 1- Note 
that the long-dashed curve shows the contribution to the B-mode from only m = lorm = — 1 (not 
the sum of m = 1 and m = —1), and similarly the short-dashed curve shows the contribution to 
the E-mode from only m = 2 or m = — 2. We find that the dominant contribution to the E-mode 
is from m = and rn = ±2, while for the B-modc it is from m = ±1. There is no B-modc for 
m = and very small E- and B-modes for m = ±1 and m = ±2, respectively. The total power 
spectrum of the E-mode, plotted as the solid curve, and calculated by summing the m = 0, ±1 and 
±2 contributions, is the same as the total power spectrum of the B-mode. The power spectra are 
multiplied by T^mb- 
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Fig. 4. — Power spectra of the second order effect for our models (note that the curve for model 
D is hidden under that for model E). For comparison, the toy model of Gruzinov & Hu (1998), 
in which the ionized regions are randomly distributed, is also shown. The heavy solid line shows 
the first-order anisotropy for J7b = 0.02 and /esc = 1-0. Note that only Cei contributes to the 
first-order effect, while in the second-order effect Cei and Cbi are equal. 
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Fig. 5. — Effect on the second-order anisotropy of varying the assumed geometry of the re-ionized 
regions. The models are all computed for l^b = 0.02 and /esc = 1-0. The line for model A is 
repeated from Fig. 4. The line labelled "uniform Xg" is for the case where the ionized fraction 
is uniform, and equal to the mean value in model A. For the line labelled "random", the total 
gas density is forced to be uniform, and spheres are put down at random positions but with the 
distribution of radii given by the galaxy formation model. The case "random (equal radii)" is the 
same, except that all spheres have the same comoving radius (0.62/i~^Mpc). Finally, for the case 
"clustered (equal radii)", variations in the total gas density are included, but all spheres have equal 
radii at a given redshift, and are placed on random dark halos. See the text for more details. 
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Fig. 6. — Effect on the secondary CMB polarization anisotropy spectrum of varying fi]-,. Tlie 
curves shown are aU computed using A-Iodel A. The sohd hue shows il^ = 0.04 while the dotted Unc 
shows ilb = 0.02. The heavy sohd hne shows the first-order anisotropy for the model Ob = 0.04 
and /esc = 1-0. The power spectrum is multiplied by Tcmb- 
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Fig. 7. — Effect on the secondary CMB polarization anisotropy spectrum of varying the escape 
fraction /esc- The curves shown are all computed using Model A and 17^ = 0.02. Escape fractions 
of 1.0,0.1 and 0.05 are shown as indicated in the figure. The heavy solid line shows the primary 
anisotropy for the model Jib = 0.02. The power spectrum is multiplied by ^cf^g. 



